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ABSTRACT 

In  studying  light  and  image  transfer  in  sea  water  the  influence  of  Fresnel  surface 
reflection  is  as  significant  as  scattering  and  absorption  phenomena.  In  these  cases  a 
knowledge  of  the  reflective  properties  of  sea  surface  at  different  wind  speeds  is  very 
important.  At  present,  little  is  published  about  these  properties.  We  present  here  results 
of  numerical  modeling  of  Fresnel  light  reflection  coefficient  of  sea  water  as  a  function  of 
solar  zenith  angle  and  wind  speed.  The  ray-tracing  computer  model  was  developed  to 
generate  wave  slopes  and  elevations.  In  order  to  generate  a  realistic  sea  surface  the  model 
used  Paul  Hwang  wave  height  spectrums.  The  final  result  of  this  paper  is  a  simple 
equation  and  very  fast  FORTRAN  code  to  calculate  Fresnel  reflection  coefficient  of  wavy 
water  surface  and  specular  part  of  remote  sensing  reflectance. 


INTRODUCTION 

The  optical  remote  sensing  reflectance  of  the  ocean  consists  of  two  parts:  one  that  is  connected 
to  the  specular  reflection  of  light  from  wavy  water  surface,  and  the  second  one  that  is  due  to  the 
diffuse  light  ascending  from  the  water.  The  information  about  the  optical  state  of  water  body  is 
contained  in  the  internal  light  coming  from  the  sea.  In  order  to  precisely  estimate  the  second 
internal  component  of  the  remote  sensing  reflection  we  need  an  approach  to  precisely  estimate  the 
first,  specular  reflected  component  of  light. 

This  presentation  explains  a  new  and  very  fast  algorithm  to  compute  specular  reflection  of  solar 
light  from  the  wavy  water  surface.  The  specular  light  reflected  from  the  water  surface  consists  of  two 
parts.  The  first  part  is  due  to  the  Fresnel  reflection  from  a  wind-roughened  surface,  and  the  second 
one  is  due  to  the  reflection  by  whitecaps.  The  Fresnel  part  of  reflection  was  computed  using  a  Monte 
Carlo  method.  The  roughness  of  water  surface  was  realistically  modeled  with  the  use  of  surface  wave 
energy  distributions  by  Paul  Hwang.  The  resulting  reflection  coefficient  of  light  by  a  water  surface 
was  obtained  by  averaging  over  billion  realizations  of  reflection  coefficient.  Results  of  these 
calculations  are  expressed  as  a  very  precise  regression  that  connects  Fresnel  reflection  coefficient  to 
the  solar  angle  and  wind  speed. 

The  proposed  algorithm  is  implemented  as  a  very  fast  FORTRAN  code  (See  Appendices  A  and  B) 
to  compute  a  specular  reflection  component  of  the  remote  sensing  reflectance. 


Presented  to  the  Seventh  International  Conference:  Remote  Sensing  for  Marine  and  Coastal 
Environments,  Miami,  Florida,  USA,  20-22  May  2002. 
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SPECULAR  PART  OF  REMOTE  SENSING  REFLECTANCE 

The  specular  part  of  remote  sensing  reflectance  consists  of  two  parts,  the  first  one, 
connected  with  Fresnel  reflection  coefficient  of  light  from  rough  (wavy)  water  surface  and  the 
second  one,  related  to  reflection  of  light  from  sea  foam  or  whitecaps  It  can  be  expressed  as  a 
following  expression: 


^5  “  ^F’  (1) 

The  reflection  by  sea  foam  can  be  regarded  as  Lambertian  and  is  expressed  as  follows: 

rsF^f(ti)Ap/7:,  (2) 

where  is  an  albedo  of  the  foam  (A^  =  0.6),  and  /(u)  is  a  portion  of  the  sea  surface  covered  with 
whitecaps  (foam),  w  is  a  wind  speed.  The  function  /(u)  is  defined  by  the  following  regression 
(Frouin,  1996): 


/(«)  = 


0.000012  u  <  9m/s, 
(0.225  w- 0.99)  u>  9m /s. 


(3) 


The  Fresnel  portion  of  the  remote  sensing  reflectance  is  defined  by  the  following  equation: 


fp  =  [i-  , 


(4) 


where  is  an  angle  of  incidence  equal  to  a  solar  zenith  angle  for  flat  water  surface.  When  wind 
speed  is  equal  to  zero,  there  are  no  foam  on  a  sea  surface  (/(w)  =  0)  and  foam-related  part  of  remote 
sensing  reflectance  is  equal  to  zero  rgp  =  0,  the  surface  is  flat  and  rp  =  Rp(Zg),  where  Rp(Zg)  is 
Fresnel  specular  reflection  coefficient  from  flat  water  surface,  defined  by  the  following  equation: 


R^p(Zs)  = 


cosZo 


sin  Ze 


cosZ^  -sin^Z, 


'5  7 


^  nl  cos  Zs-^jnl-  sin^  Z^  1 
,  nl  cos  Zs  +  ^Inl-sin^Zs  j 


(5) 


here  -1.341  is  a  refractive  index  of  water. 

When  wind  started  to  blow,  it  roughens  the  water  surface  and  the  situation  becomes  more 
complex:  we  need  to  average  expression  (5)  over  incidence  angles.  There  are  different  approaches  to 
accomplish  this  procedure:  1)  we  can  use  a  Cox  and  Munk  (1954)  sea  water  slope  distribution,  or  2) 
we  can  use  Monte  Carlo  method  to  generate  wave  slopes  using  wave  energy  spectrums  taken  from 
experimental  measurements  (Haltrin,  McBride  III,  and  Weidemann,  2000).  The  second  approach  is 
more  complex  and  time-consuming,  but  it  gives  better  results,  especially  for  remote  sensing 
applications.  To  compute  specular  part  of  remote  sensing  reflectance  in  this  paper  we  used  the 
second  approach. 
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COMPUTATIONS  OF  FRESNEL  REFLECTION  COEFFICIENT 

To  generate  sea  water  surface  slopes  we  used  a  Monte  Carlo  program  described  in  Haltrin, 
McBride  III,  and  Arnone  (2001).  The  values  of  sea  surface  orientation  allow  us  to  compute  light 
reflection  coefficient  in  each  pixel  of  the  sample  surface.  Averaging  over  MxM  spatial  and  L 
temporal  realizations  gives  as  an  average  Fresnel  reflection  coefficient  for  each  value  of  wind  speed 
u  and  zenith  angle  Z^.  To  compute  actual  values  of  Fresnel  reflection  coefficient  we  used  the  Paul 
Hwang  (1997)  spectrum.  This  spectrum  is  specifically  tailored  for  remote  sensing  problems  to 
produce  correct  values  of  mean  square  slopes  of  ocean  waves. 

The  generated  realizations  of  Fresnel  reflection  coefficients  have  been  averaged  over 
100x100  pixels  of  sample  sea  surface  areas  and  80  time  realizations  to  produce  resulting  angular 
distributions  of  Fresnel  reflection  coefficient.  The  values  of  elevations  and  orientations  in  each  pixel 
have  been  obtained  using  one  million  computations,  1000  realizations  for  flat  sine  waves  and  1000 
significant  points  of  wave  number  k  in  energy  spectrum  range.  To  generate  random  numbers  we  used 
Mersenne  Twister  random  number  generator  (Matsumoto  and  Kurita,  1992,  1994)  capable  to 
produce  evenly  distributed  random  numbers  in  a  cube  of  626  dimensions.  So  each  value  of  Fresnel 
reflection  coefficient  for  any  value  of  wind  speed  and  zenith  angle  represents  an  average  value  of  800 
billion  individual  computations. 


FRESNEL  REFLECTION  COEFFICIENT  BY  WAVY  WATER  SURFACE 

The  results  of  extensive  Monte  Carlo  computations  are  compactly  represented  as  a  following 
regression  equation  (Haltrin,  McBride  III,  and  Amone,  2001): 

R,{u,Zs)  =  «o(“)  +  +  Rl(Zs)[a,{u)  +  fl3(M)J?;(Z,)]},  (6) 

where  /?^(Z^)  is  a  Fresnel  reflection  coefficient  of  flat  water  surface  given  by  Eq.  (5),  and  wind- 
speed-dependent  coefficients  a.{u)  are  given  by  the  following  equations: 

a^{u)  =  0.001(6.944831- 1.912076m +0.03654833m^),  =  0.9997,  (7) 

a,{u)  =  0.7431368 +0.0679787  M- 0.0007171  M^  r'  =  0.9996,  (8) 

^^(m)  =  0.5650262 +0.0061502M-0.0239810M^+0.0010695M^  =  0.9995,  (9) 

a^{u)  =  -0.4128083-0.1271037m+0.0283907m"-0.0011706m\  r'  =  0.9991.  (10) 


Equations  (1)“(6)  represent  a  complete  algorithm  to  compute  the  specular  part  of  remote  sensing 
reflectance.  This  algorithm  is  programmed  in  FORTRAN  and  presented  in  Appendices  A  and  B. 
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APPENDIX  A;  PROGRAM  RSSPEC.F 


I 

! 

I 

! 

! 


program  rsspec 


Vladimir  I.  Haltrin  <haltrin@nrlssc.navy.mil> 
last  modification:  February  20,  2002 
********************************************************************** 


vwk 
vwm 
f,  ff 

Zs ( 1 :nSun) 
nWat 


=  wind  speed  over  the  sea  surface  in  knots 

=  wind  speed  over  the  sea  surface  in  m/sec 

=  fraction  of  sea  surface  covered  by  foam 

=  a  massive  of  solar  zenith  angles 
-  refractive  index  of  sea  water 
********************************************************* 


implicit 

integer 

parameter 

real 

real 

real 

logical 


none 

nSun,Ns,Nw,Nz, j ,k 
(Nw=20,Nz=90) 

nWat ,R, Afoam,Pi,rd,ws ,dt , f f ,phi, frO^rsa, f (Nw) 
windspeed(Nw) ,vwk(Nw) ,vwm{Nw) , Zs{Nz) ,RfrFl(Nz) , Fresnel 
RfrWn(Nw,Nz ) ,Rrs (Nw,Nz ) ,Ff oam, FresRwind 
knots 


open (11, f ile= ’ rsspec .in', status= ' old ' ) 
read (11,*)  Afoam 
read(ll,*)  nWat 
read (11,*)  knots 
read(ll,*)  Ns 

read(ll,*)  (windspeed( j ) , j=l,Ns) 
read (11,*)  nSun 

read(ll,*)  (Zs(k),  k=l,  nSun)  !  solar  zenith  angles  in  degrees 
close ( 11 ) 


Pi  =  4.*ATAN(1. ) 
rd  =  Pi/180. 


if  (knots)  then 
do  j=l,Ns 
vwm(  j )  - 
vwk(j)  = 
end  do 
else 

do  j=l,Ns 
vwm(  j  )  “ 
vwk ( j )  = 
end  do 
end  if 


0 .515*windspeed( j ) 
windspeed( j ) 


windspeed( j ) 
windspeed( j ) /O .515 


do  k  =  1,  nSun  !  solar  zenith  angle  loop 

phi  =  rd*Zs(k) 

frO  =  Fresnel (nWat, phi)  !  Fresnel  reflection  coefficient 

RfrFl(k)  =  frO 

do  j=l,  Ns  !  windspeed  loop 

ws  =  vwm( j ) 
ff  =  Ffoam(ws) 
f(j)  =  ff 

R  =  FresRwind( f rO ,ws ) 

RfrWn(j,k)  =  R 

Rrs(j,k)  =  ff*Afoam/Pi  +  (l.-ff)*R 
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end  do 
end  do 

call  Rout( 'RfrWnd.out' ,Afoam,nWat,Nw,N2,Ns,nSun 
&  ,f ,vwk,vwm/Zs,RfrFl,RfrWn) 

call  Rout( 'RSspec.out' , Afoam,nWat,Nw,Nz,Ns,nSun 
&  , f ^ vwk^vwm, Zs^RfrFl/Rrs ) 

end 

I  ********************************************************************** 

real  function  Fresnel (nWat/ phi) 


i  phi  =  angle  (in  radians)  at  which  light  is  incident  on  sea 

!  aRef  =  angle  of  refraction  (in  radians) 

!  Rpar  =  Fresnel  refrlection  coefficient  for  parallel  polarization 

!  Rper  =  Fresnel  refrlection  coefficient  for  perpend,  polarization 

I  ********************************************************************** 

implicit  none 

real  nWat , phi ,  aRef , aDif , aSum, Rpar , Rper 

if  (phi  .ne.  0.)  then 

aRef  =  ASIN(SIN(phi)/nWat) 

aDif  =  phi-aRef 

aSum  =  phi+aRef 

Rpar  =  TAN(aDif )/TAN(aSum) 

Rper  =  SIN(aDif )/SIN(aSum) 

Fresnel  =  0 . 5* (Rpar*Rpar+Rper*Rper ) 
else 

aSum  =  (nWat~l . ) / (nWat+1 . ) 

Fresnel  =  aSum*aSum 
end  if 

return 

end 

I  ********************************************************************** 

real  function  FresRwind( frO ,ws ) 


Calculates  Fresnel  reflection  coefficient  of  wavy  water  surface 
for  the  case  of  Paul  Hwang  wave  energy  spectrum  distribution 
frO  =  Fresnel  reflection  coefficient  of  a  flat  water  surface, 
ws  =  windspeed  in  m/sec,  0  <=  ws  <=  12  m/s. 

See;  V.  I.  Haltrin,  W.  E.  McBride  III,  and  R.  A.  Arnone,  "Spectral 
approach  to  calculate  specular  reflection  of  light  from  wavy  water 
surface,"  -  pp.  133-138  in  Proceedings  of  D.  S.  Rozhdestvensky 
Optical  Society:  International  Conference  Current  Problems  in  Optics 
of  Natural  Waters  (ONW'2001),  St.  Petersburg,  Russia,  2001. 
********************************************************************** 

implicit  none 

real  fr0,ws,  a0,al,a2,a3 

aO  =  0.001*(6.944831+ws*(-l. 912076+0. 03654833*ws)) 
al  =  0.7431368+ws*(0. 0679787-0. 0007171*ws) 

a2  =  0.5650262+ws*(0.0061502+ws*{-0. 023981+0. 0010695*ws)) 
a3  =  -0. 4128083+ws*(-0.1271037+ws*(0. 0283907-0. 0011706*ws) ) 

FresRwind  =  a0+fr0*(al+fr0*(a2+a3*fr0) ) 

return 
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end 


I  ********************************************************************** 

real  function  Ffoam(ws) 

I  - 

!  Calculates  a  fraction  of  water  surface  covered  by  foam 

I  ws  =  windspeed  in  m/sec,  0  <=  ws  <=  12  m/s. 

!  R.  Frouin,  Ocean  Optics  XIII,  Halifax,  Canada,  unpublished. 

I  ********************************************************************** 

implicit  none 
real  ws,f 

f  =  ws*ws*ws 

f  =  1.2E-5*f*{ws**0.3) 

if  (ws  .gt.  9.)  f  -  f * ( 0 . 225*ws-0 - 99 ) 

if  ( f  .gt .  1 . )  f  =  1 . 

Ffoam  =  f 

return 

end 


I 


subroutine  Rout ( f name, Afoam,nWat,Nw,Nz , Ns ,nSun 
&  , f ,vwk, vwm, Zs,RfrFl,Rf ) 

********************************************************************** 


implicit 

integer 

real 

real 

character 


none 

Nw,Nz,Ns,nSun,k, j 
Afoam,nWat,f (Nw) ,vwk(Nw) ,vwm(Nw) 
Zs(Nz) ,RfrFl(Nz) ,Rf (NW,NZ) 
tb, fname*10 


tb  =  CHAR(9) 

open (21,  file=fname,  status='new' ,recl=250) 
write (2 1,44)  Afoam 
write (2 1,55)  nWat 

write(21,66)  '  f : ’ , tb, 0 . , (tb,  f(j),j=l,Ns) 

write (2 1,66)  ' ZsNvwk: ' ,tb,0. , (tb,vwk( j ) , j=l,Ns) 
write (2 1,66)  ' - 


write (2 1,66)  'ZsXvwm: ’ ,tb,0., (tb, vwm( j ) , j^l ,Ns ) 
do  k  =  1,  nSun 

write (2 1,77)  Zs (k) , tb,RfrFl(k) , (tb,Rf ( j ,k) , j=l ,Ns ) 
end  do 
close ( 21 ) 


44  format (X, ' Foam  Albedo  =  ',f7.5) 
55  format ( X, 'Water  refraction  coefficient  =  ',f7.5) 
66  format (a7,  21 (al, f 10.5 ) ) 

77  format(f6.3,21(al,fl0.7) ) 


return 

end 
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APPENDIX  B:  INPUT  FILE  RSSPEC.IN 


0.60 

< — 

Afoam 

,  foam 

1.341 

< — 

nWat , 

water 

. false. 

< — 

knots 

(if  .1 

6 

< — 

Ns ;  ^ 

2. 

4.00 

6. 

8. 

10. ( 

22 

< — 

nSun; 

o 

o 

05. 

10. 

15 

20. 

50, 

55. 

60. 

65. 

70. 

12. 

sun  zenith  angles  in  degrees,  ang(l;nSun): 
25.  30.  35.  40.  45. 

75.  77.5  80.  82.5  85.  87.5  89.9 
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